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Abstract 



Let H be the space of all Hermitian matrices of infinite order and U{oo) be the 
inductive limit of the chain C C/(2) C . . . of compact unitary groups. The 

group U{oo) operates on the space H by conjugations, and our aim is to classify 
the ergodic U (oo)-invariant probability measures on H by making use of a general 
asymptotic approach proposed in Vershik's note [V]. The problem is reduced to 
studying the limit behavior of orbital integrals of the form 



where A is a fixed cxo x oo Hermitian matrix with finitely many nonzero entries, 
fin is a [/ (n) -orbit in the space of n x n Hermitian matrices, Mn is the normalized 
?7(n)-invariant measure on the orbit Qn, and n — > oo. 

We also present a detailed proof of an ergodic theorem for inductive limits of 
compact groups that has been announced in [V]. 

There is a remarkable link between our subject and Schoenberg's [S2] theory of 
totally positive functions, and our approach leads to a new proof of Schoenberg's 
[S2] main theorem, originally proved by function-theoretic methods. 

On the other hand, our results have a representation-theoretic interpretation, 
because the ergodic t/(oo) -invariant measures on H determine irreducible unitary 
spherical representations of an infinite-dimensional Cartan motion group. 

The present paper is closely connected with a series of articles by S. V. Kerov 
and the authors on the asymptotic representation theory of "big" groups, but it 
can be read independently. 



Both authors were partiahy supported by the International Science Foundation under Grant 
MQVOOO. The first-named author was also partially supported by the Russian Foundation for 
Basic Research under Grant 95-01-00814. 
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To the memory of F. A. Berezin 



§0. Introduction 



1. Ergodic measures. The description of ergodic invariant measures for group 
actions is a traditional problem of ergodic theory. It is well known that it is not 
always that this problem can be solved in a satisfactory way. It is a surprising 
fact that for certain infinite-dimensional (or "big"^) groups, there exist nice actions 
whose ergodic measures can be completely described. The simplest examples are 
given by the classical theorems due to B. de Finetti and I. J. Schoenberg. 

We recall that, by de Finetti's theorem, the ergodic measures on a product 
space X°°, invariant under the group of permutations of the copies of X (the 
infinite symmetric group), are exactly the product measures with identical factors. 
Schoenberg's theorem [SI] states that the ergodic 0(oo)-invariant measures on the 
space M°° are exactly the Gaussian product measures (by 0(oo) we denote the 
inductive limit group linjO(n)).^ 

These classical examples and a number of more complicated ones were discussed 
in the note [V] by one of the authors. In that note, a general "ergodic method" 
for inductive limits of compact groups was proposed. This method was further 
developed in the papers [VKl, VK2, KV] by Vcrshik and Kcrov. The aim of the 
present paper is to give a detailed exposition of the method on a model example 
considered in [V]. 

Let H{n) denote the space of n x n complex Hermitian matrices, and let H = 
^m.H{n) be the space of all infinite Hermitian matrices. Let U{oo) = lin^ U (n) 
be the group of infinite unitary matrices u = [uij] such that Uij = 5ij when i + j 
is large enough. The group U{oo) operates on the space H by conjugations, and 
we are interested in the class M. of all ergodic U ((X))-invariant Borel probability 
measures on H. 

Further, let H{oo) = linji?(n) be the space of oo x oo Hermitian matrices with 
finitely many nonzero entries. The spaces H{oo) and H are in a natural duality, 
and any measure on H is uniquely determined by its characteristic function (Fourier 
transform), which is a function on H{oo). 

Classification Theorem. The characteristic functions of the measures M ^ M. 
are exactly those of the form 



where 71,72,2:1,2:25 ■ ■ ■ , the parameters of the measure M, are real numbers such 



The Classification Theorem implies that any measure M G can be written as 
the convolution product of a Gaussian measure and a (finite or countable) family of 

■""The term "big groups" was suggested by one of us in [V]. This term has no rigorous definition, 
it can mean "infinite-dimensional" or "out of the class of locally compact groups" or something like 
that. For instance, the group S{oo) := finj 5(n), the inductive limit of finite symmetric groups, 
is a discrete, hence a locally compact group, but its properties are very similar to that of the 
infinite-dimensional group U{oo) :— lirt^ U(n), so we prefer to rank 5(oo) among the "big" groups. 

^ About this theorem, see also Berg-Christensen-Ressel [BCR]. 
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non-Gaussian "elementary" ergodic measures. The latter are essentially supported 
by rank one matrices and are related to Wishart distributions, well-known in mul- 
tivariate statistical analysis. Let us emphasize that in this theorem, the structure 
of the answer turns out to be much more complicated than in de Finetti's and 
Schoenberg's [SI] theorems: up to trivial exceptions, the ergodic measures M e A4 
are neither product nor Gaussian measures. 

Note that formula (0.1) can be rewritten in the following form: 

f{A)= n ^(«)' AeH{oo), (0.2) 

aeSpec(A) 

where Spec (A) is the collection of eigenvalues of A (taken with their multiplicities) 
and 

F(a) = e»^i«-T'2'*'/2 rr , aeR. (0.3) 

j^J-^ 1 - ixka 

The function F{a) has a simple meaning. Let us regard the matrices B E H as 
random matrix variables defined on the probability space (if, M). Then F{a) is the 
characteristic function for the distribution of any diagonal entry Sjj, i = 1, 2, . . . 

As shown by Pickrell [Pi2], the Classification Theorem can be derived from a 
deep function-theoretic result due to Schoenberg [S2]. In the present paper, we 
prove the Classification Theorem in an entirely different way: our method consists 
in studying how an ergodic measure M G is approximated by finite-dimensional 
"orbital measures" . (By an orbital measure we mean the U (n)-invariant probability 
measure supported by a t/(n)-orbit Q.^, C H{n), where n = 1,2,....) 

We would like to emphasize that such an approach not only leads to the descrip- 
tion (0.1) of ergodic measures, but also provides us with additional information 
about them. Namely, we find a characterization of those sequences {M^} of orbital 
measures that weakly converge, as n — > oo, to an ergodic measure M E M., so we 
can understand how M "grows up" from orbital measures. In particular, we can re- 
late the parameters 71, 72, iCi, iC2) • • • to the asymptotic behavior of the eigenvalues 
of large Hermitian matrices. 

A characterization of the convergent sequences is the main result of the 
paper — Theorem 4.1. The statement is too long to be reproduced here in detail, 
but roughly speaking, the picture is as follows. Given n = 1,2,..., pick any matrix 
from the orbit Qn supporting Mn and represent the collection of the eigenvalues of 
that matrix as an n-point configuration on the real line. Then this configuration, 
being contracted with scaling factor 1/n, must converge, as n — » 00, to a (countable) 
point configuration {xi, X2, ■ ■ ■} C M. (Here the x^s coincide with the parameters 
appearing in (0.1); in fact, there are also additional restrictions on the growth of 
the orbits related to the parameters 71 and 72.) 

2. Relationships with representations and totally positive functions. We 

are interested in the measures of class for several reasons: 

1) (The initial motivation of [V].) The action of the group L'^(oo) on the space 
if is a natural example of a "big" group action. As compared with the action of 



^This implies that the whole distribution of the matrix entries is uniquely determined by the 
distribution of any diagonal entry! Note that this property holds for ergodic measures M only: it 
follows from the so-called multiplicativity theorem, see Theorem 2.1 below. 
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the group 0(oo) on the space M°° (Schoenberg's case [SI]) or the action of U{oo) 
on C°°, this example is of the next level of complexity. So, from the viewpoint of 
ergodic theory, it is interesting to find ergodic measures and to compare the result 
with the one in Schoenberg's case. 

2) The ergodic measures of class A4 determine some irreducible unitary repre- 
sentations of a certain infinite-dimensional group G{oo). 

3) These measures are closely related to Schoenberg's totally positive functions. 

4) They are also related to harmonic analysis on the group U{oo) (cf. [KOV]). 
We shall briefly comment on items 2) and 3). 

The group G{oo) is a model example of an infinite-dimensional Cartan motion 
group. It can be defined as the inductive limit G{oo) = linjG(n), where G{n) stands 
for the scmidircct product U{n) tx H{n). There is a natural bijective correspondence 
between ergodic t/(oo)-invariant probability measures on H and irreducible unitary 
representations of the group G{oo) that are spherical with respect to the subgroup 
U{oo) C G(oo). Under this correspondence, the spherical functions of G(oo), when 
restricted to H(oo) := linjiy(n), coincide with the characteristic functions of the 
ergodic measures. 

It should be noted that all the results and constructions concerning the ergodic 
measures of class /A can be translated into the language of representation theory. 
In particular, weak convergence of measures turns into uniform compact conver- 
gence of spherical functions, and the problem that is solved in the present paper is 
a particular case of the following general problem in representation theory of "big" 
(inductive limit) groups: Given a (irreducible) unitary representation of an induc- 
tive limit group Q — lin^^(n), study its approximation by (irreducible) unitary 
representations of the growing subgroups Q{n) as n — > oo. 

(About various asymptotic results in representation theory of "big" groups, see 
the scries of papers by the authors and S. V. Kerov, [VK1-VK3, Ke, KV, 01-07].) 

A remarkable feature of representation theory of "big" groups is its connection 
with analytic problems of total positivity theory (see Thoma [T], Boyer [Bo], Ver- 
shik and Kerov [VKl, VK2], Voiculescu [Vol]; for a systematic exposition of total 
positivity theory, see Karlin's fundamental monograph [K]). 

To describe this connection, we need the definition of totally positive (TP) func- 
tions: these are real- valued functions (f{t) on M such that for any n and any choice 
of real numbers ti < ■ ■ ■ < tn and si < • ■ • < the determinant of the n x n ma- 
trix [(f{ti — Sj)] is nonnegative. Similarly, one also defines two-sided and one-sided 
TP sequences: these are TP functions defined on the 1-dimensional lattice Z C IR 
or its nonnegative part respectively. 

It turns out that the Fourier transform of TP functions on M leads exactly 
to functions F of the form (0.3): this claim is Schoenberg's main theorem in [S2]. 
Thus, there is a correspondence M <-> between ergodic measures M e A4 and 
TP functions ip (normalized by the condition J ip{t) dt — 1) that can be stated as 
follows: (f{t) dt coincides with the distribution of (any) diagonal entry Ba, where 
B stands for a random Hermitian matrix distributed according to M. 

Further, there exists a similar correspondence between characters of the group 
U{oo) or of the infinite symmetric group S'(oo) and two-sided or one-sided TP 
sequences, respectively. 

This correspondence can be used in two directions: 

On the one hand, old theorems in total positivity theory, obtained by ana- 
lytic tools in [ASW, El, E2, S2, T, K], can be applied to representation-theoretic 
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problems (classification of spherical functions or characters). Such an approach is 
adopted by Thoma [T], Boyer [Bo], and PickrcU [Pi2]. 

But on the other hand, if we can classify spherical functions or characters inde- 
pendently, then we can prove some theorems on TP functions or TP sequences in 
a new way. 

It is the second approach that is adopted in Vershik and Kerov's papers [VKl, 
VK2], and also in the present paper: as a corollary of our main result, we obtain a 
new derivation of Schocnbcrg's classification for TP functions on the real line. 

Also note a recent paper by Okounkov [Ok] , where a different (direct representa- 
tion-theoretic) method is used; the result of [Ok] gives yet another way to classify 
one-sided TP sequences. 

3. The method. Now let us describe our techniques in more detail. 

The starting point of our approach is a general approximation theorem for er- 
godic measures, see Theorem 3.2 below^. When applied to our concrete situation, 
it implies that any ergodic measure M E M. can be approximated by a sequence 
{Mn} of orbital measures (see Theorem 3.3 below). Then the main problem is to 
understand what sequences {M^} are weakly convergent, as n — > oo, and what are 
their limits. 

An attempt to do this was made in [V], but the solution proposed there was 
incomplete, because of a gap in the calculations. Nevertheless, the method itself 
was correct, and a refinement of the calculations leads to the right result. 

To study the weak convergence of orbital measures M„, we deal with their char- 
acteristic functions /„, 



where Qn C H{n) is the t/(n)-orbit carrying Mn. We calculate the Taylor decompo- 
sition at the origin of a characteristic function /„ and then analyze the asymptotics 
of its Taylor coefficients as n — > oo. Note that these Taylor coefficients are nothing 
but the moments of the measure M„. (In fact, due to the symmetry of fn, it is 
convenient to rewrite its Taylor decomposition as a series of Schur polynomials.) A 
non evident fact is that weak convergence of orbital measures can always be con- 
trolled by the moments. Such a phenomenon was first discovered in [VK2]; there 
it was applied to an allied problem: classifying characters of the group U{oo). 

Note that we are able to generalize our results to the spaces of real symmetric 
and quaternionic Hermitian matrices. Instead of Schur polynomials, we must then 
use Jack symmetric polynomials. As a further generalization, one could consider 
all matrix spaces of the form H = ^miy(n), where H{n) ranges over one of the 
10 series of classical symmetric spaces of Euclidean type, the parameter n being 
the rank of the space (see [05, Pil] and especially [Pi2] for a discussion of these 
spaces H). Characteristic functions of orbital measures on the spaces H{n) are 
sometimes called generalized Bessel functions; for 3 of the 10 spaces H, they can 
be expressed in terms of elementary functions, and for the 7 other spaces, these are 
certain multidimensional special functions. The problem consists in studying their 
limiting behavior as n — > oo. We conjecture that our approach can be transferred 
to all the spaces H. 




•71 



(0.4) 



'^There is also another version of this theorem that deals with spherical functions, see Theorem 
3.5 below. 
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4. Contents. The present paper is organized as follows. 

In §1, we introduce ergodic measures and explain their relationship with spherical 
unitary representations. 

In §2, we formulate the so-called Multiplicativity Theorem for characteristic 
functions of ergodic measures (Theorem 2.1). This important result states that an 
U (oo)-invariant probability measure M on H is ergodic if and only if its character- 
istic function is multiplicative in a certain sense; it follows that the set of ergodic 
measures is closed under convolution. Then we state the classification result (Theo- 
rem 2.9), which shows that any ergodic measure is a convolution product of certain 
"elementary" measures. We also discuss a number of corollaries of these theorems. 

Section 3 is devoted to a proof of a general approximation theorem that has been 
announced in [V] (Theorem 3.2 below). This result may be viewed as an ergodic 
theorem for actions of general inductive limits of compact groups. 

In §4, we state the main result of the paper, Theorem 4.1, and outline its proof. 

Section 5 contains a preliminary result, needed for the proof of Theorem 4.1. 
There we decompose the characteristic function of an orbital measure into a series 
of Schur symmetric functions.^ 

In §6, we prove Theorem 4.1. 

In §7, following Pickrell's arguments in [Pi2], we establish the equivalence of two 
classification problems: that of ergodic invariant measures M G J\4 and that of 
totally positive functions (p on R. Note that our own contribution to the results of 
this section is very modest. Here we only aimed to clarify some technical details of 
the correspondence M <-> <^ and to explain how our main result implies Schoenberg's 
classification.^ 

In the final §8, we explain a connection between our main theorem and the main 
result of the remarkable work [CS] by Curry and Schoenberg. 

A preliminary version of the present paper appeared as preprint [OV]. 

5. Acknowledgements. Already in the sixties, the second named author (A. V.) 
had discussions with F. A. Berezin about the subject of the present work and its 
links with mathematical physics. 

The preparation of the preliminary version [OV] of the paper was completed while 
the first named author (G. O.) was a guest of the Department of Mathematical 
Sciences, University of Tokyo. It is a pleasure to thank Masatoshi Noumi and 
the Department of Mathematical Sciences for the kind invitation to Tokyo and 
hospitality. 

Special thanks are due to Andrei Okounkov who brought to our attention the 
important relation with the work [CS] by Curry and Schoenberg. 

The authors also acknowledge the financial support by the International Science 
Foundation and by the Russian Foundation for Basic Research during the final 
phase of our work. 

§1. Ergodic measures and spherical representations 
Let H{n) denote the real vector space formed by complex Hermitian n x n 



°It would be interesting to find analogues of this result for all generalized Bessel functions. 

^Note, however, that we are not completely satisfied by the way we obtain Schoenberg's the- 
orem. We think there should exist a more conceptual derivation of this result by approximation 
methods. A similar remark can be made d propos TP sequences as well. 
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matrices, n = 1,2,.... There is a natural embedding 



H{n) ^ H{n+1), 

and we denote by H{oo) the corresponding inductive hmit space linji7(?7,). Then 
H{oo) is identified with the space of infinite Hermitian matrices with a finite number 
of nonzero entries. We equip if (oo) with the inductive limit topology. In particular, 
a function /: H{oo) — > C is continuous if its restriction to H{n) is continuous for 
any n. 

Let H stand for the space of all infinite Hermitian matrices. For A E H and 
n = 1, 2, . . . , we denote by On{A) G H{n) the upper left n x n corner of A. Using 
the projections 6^: H ^ H{n), n = 1, 2, . . . , wc may identify H with the projective 
limit space ^mH{n). We equip H with the corresponding projective limit topology. 
In particular, H is a Borel space. 

There is a natural pairing 



A 




H{oo)xH^R, {A, B) ^ tr{AB) . 



DCX) 



Using it, we may regard H as the algebraic dual space of H{oo). 

Note that H can be identified, in an obvious manner, with the space 
M X R X ■ ■ ■ . Under this identification, H{oo) C H turns into Mg" •= Un>i '^^^ 
the pairing defined above becomes the standard pairing between and R°°. 

Given a Borel probability measure M on H, we define its Fourier transform, or 
characteristic function, as the following function on H{(X)): 

f{A) = fM{A)= I e''<^^^M{dB). (1.1) 

JH 

We need the following statement: 

Proposition 1.1. The Fourier transform (1.1) establishes a bijective correspon- 
dence between Borel probability measures M on H and continuous positive definite 
normalized functions f on H{oo). 

Here "normalized" means that /(O) = 1. 

Proof. This is an immediate corollary of Kolmogorov's consistency theorem and 
Bochner's theorem. Indeed, let us identify {H{oo),H) with (Rg°,M°°). Then Kol- 
mogorov's theorem implies that Borel probability measures on are just the 
projective limits of probability measures on the ]R"^'s, and Bochner's theorem al- 
lows us to restate this fact in terms of characteristic functions. □ 

Let U{n) be the group of unitary n x n matrices, n = 1, 2, ... . For any n, we 

\ u 1 

embed U{n) into U{n+ 1) using the mapping u ^ ^ ^ . Let U{oo) = Iwcj. U (n) 

denote the corresponding inductive limit group. We regard U{oo) as the group of 
infinite unitary matrices u = [uij]fj^i with a finite number of entries Uij ^ Sij. The 
group U{oo) acts by conjugations both on H{oo) and H, and the pairing between 
these two spaces is clearly U (oo) -invariant. 

Let us recall some basic facts concerning ergodic measures. 



Definition 1.2. Let A" be a Borel space, let ^ be a group of Borel transformations 
of A", and let M be a ^-invariant probability Borel measure on X. 

(i) A Borel subset y G X is said to be Q-invariant mod if M{{gy)Ay) = 
for any g E Q. 

(ii) The measure M is said to be ergodic with respect to Q if the M-volume of 
any ^-invariant mod Borel subset y C X is equal either to or 1. 

Proposition 1.3. The following conditions on a Borel probability measure M on 
X are equivalent: 

(i) M is ergodic; 

(ii) M is an extreme point of the convex set formed by all Borel probability 

measures on X; 

(iii) the subspace of Q-invariant vectors in the Hilbert space L'^{X,M) is ex- 
hausted by constant functions. 

Proof. See, e.g., Phelps [Ph, Proposition 10.4]. □ 

For compact group actions on locally compact spaces, ergodic measures coincide 
with invariant measures supported by the orbits, so that for such actions, classifying 
orbits and classifying ergodic measures are equivalent problems. However, for "big" 
groups like U{oo), these two problems are quite different: the first one seems to be 
out of reach, while the second one has, in certain cases, a nice solution. According 
to a general principle of ergodic theory, ergodic measures may be viewed as true 
substitutes of orbits. 

Let Ai denote the set of all ergodic C/(oo)-invariant Borel probability measures 
M on H, and let J-' denote the set of the corresponding characteristic functions 
/ = /m- 

The following claim, which follows from Propositions 1.1 and 1.3, gives a useful 
characterization of the class J^: 

Proposition 1.4. The functions f ^ T are exactly the extreme points of the convex 
set formed by all continuous U{oo) -invariant positive definite normalized functions 
on H{oo). □ 

Let G{n) = U{n) k H{n) be the semidirect product of U{n) and the additive 
group of the vector space H{n). The elements g G G{n) are the pairs (w, A) e 
U{n) X H{n) with the multiplication rule 

(tt. A) ■ {v, B) = {uv, v~'^Av + B) . 

In the same way, we define the group G(oo) = [/(oo) k H[oo). The group G{po) 
can be also viewed as the inductive limit group linj(j(?7.), and we equip it with the 
inductive limit topology. Using the embeddings 

u^{u,Q) and A^{l,A), ueU{oo), A e H{oo), 

we may identify U{oo) and H{oo) with the corresponding subgroups in G{oo). 

Let T be a unitary representation of G{oo) in a Hilbert space TC{T) (we tacitly 
assume T is continuous with respect to the inductive limit topology of the group 
G{oo)). Then T is said to be spherical if it is irreducible and the subspace T-C{T)^^°°^ 
of C/(oo)-invariants in 7i{T) is nonzero. 

9 



Suppose T is spherical. It is then known that the space H(T)'^(°°) is one- 
dimensional (see Olshanski [05, Theorem 23.6]). A vector h E H(T)^(°°) of norm 
1 is called a spherical vector of T, and the corresponding matrix element 

<PTi9) = iTig)h,h), geGioo), 

is called the spherical function of T. Since (pT does not change when h is multiplied 
by a complex number of absolute value 1, (fiT is an invariant of T; moreover, it 
uniquely determines T. Further, as ipT is bi-invariant with respect to the subgroup 
C/(cxd) c G'(cxd), it is uniquely determined by its restriction (fT\H{oo) to H{oo), 
which is a continuous C/(oo) -invariant function. 

Proposition 1.5. There is a natural bijective correspondence T ^ M between 
{equivalence classes of) spherical representations T of the group G{oo) and ergodic 
measures M G A4. Under that correspondence, the functions ipT\H{oo) coincide 
with the characteristic functions /m- Given M E Ai, the corresponding represen- 
tation T can he realized in the Hilbert space L^{H, M) as follows: 

{T{u)-^){B) = ^{u-^Bu), u E U{oo) C G(oo), 

(T(A)*)(S) = e^^'^(^^)*(S), A E H{oo) C G'(oo), 

where G L^{H,M) and B E H. Finally, in this realization, the spherical vector 
is the constant function ^oiB) = 1. 

Proof. Indeed, it is enough to remark that the spherical functions, when restricted 
to H(oo), are just the extreme [/(cxo)-invariant continuous positive definite func- 
tions, normalized at G H{oo), and then apply Proposition 1.4. □ 

§2. The Classification Theorem 

Denote by -D(oo) the subspace of diagonal matrices in H[oo). An element of 
D{oo) will be written as diag(ai, a2, . . .), where ai, a2, • • • G R and a/. = for k 
large enough. Since any matrix in H{oo) can be diagonalized under the action of 
the group U{oo), a t/(oo)-invariant function / on H{oo) is uniquely determined by 
f\D{oo), the restriction of / to D{oo). 

Assume / is an C/(oo) -invariant function on H{oo), /(O) = 1. Let us say that / 
is multiplicative if 

/(diag(ai, a2, . . . )) = F{ai)F{a2) • • • , (2.1) 

where F is a function on R such that F{0) = 1. In other words, 

f{A)^F{ai)F{a2)-- - , AeH{oo), (2.2) 

where 01,02,... are the eigenvalues of A. Note that F{a) = /(diag(o, 0, 0, . . . )), 
i.e.,F = f\H{l). 

Theorem 2.1 (Multiplicativity Theorem). Let f be a continuous U {00) -invariant 
positive definite function on H{oo), /(O) = 1. Then f is extreme {i.e., f E if 
and only if f is multiplicative. 

Proof. See, e.g., Olshanski [05, Theorem 23.8]. There the Multiplicativity Theorem 
is stated for allied groups (like GL{n, C)) but the proof is exactly the same. □ 
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Not that there are many theorems of this type and a lot of different methods 
to prove them: see Thoma [T], Ismagilov [II] and [12] (Ismagilov's method is also 
explained in Olshanski [04, section 2.5]), Nessonov [N], Olshanski [01], Pickrell 
[Pil], Vershik-Kerov [VK4], Voiculescu [Vol] and [Vo2], Stratila-Voiculescu [SV]. 

Theorem 2.1 implies that a function / e is uniquely determined by the function 
F := f\H{l), which is a continuous function in one real variable. Note that 

F{a) = Fuia) = I e'"^'' M{dB), aeR, BeH, (2.3) 
Jh 

where M e M corresponds to /. 

Definition 2.2. Let J-'i — {F} denote the class of all functions on R of the form 
F = f\H{l) where / ranges over JF. In other words, a function F on M belongs to 
J^i if and only if F is continuous, F{0) = 1, and the corresponding function (2.2) 
on H{oo) is positive definite. 

Clearly, the classification of the measures M e 7W is reduced to that of the 
functions F e J^i. 

Theorem 2.1 has a number of corollaries. 

Corollary 2.3. The class J-'i is stable under pointwise multiplication. □ 

This implies that the class T is also stable under multiplication, and the class 
A4 is stable under convolution. 

Corollary 2.4. If a sequence Fi, F2, • ■ ■ & pointwise converges to a continuous 
function F on W, then F & !Fi. □ 

Corollary 2.5. For a real 7, the Dirac measure concentrated at 'j • 1 & H belongs 
to Ai . The corresponding function Fe.Fi is 

F{a) = e*^", aeR. (2.4) 

Proof. Since scalar matrices are C/(cxd) -invariant, our measure is invariant. It is 
clearly ergodic. It is evident that the corresponding function from Fi is given by 
(2.4). □ 

Corollary 2.6. Given 7^0, let M be the Gaussian distribution on H such that, 
for a matrix B & H, the diagonal entries B^ and the off-diagonal entries Re By, 
Im Bij, i < j, are independent Gaussian variables with mean and variance 7. 
Then M e A^, and the corresponding function F e Fi is 

F(a) = e-^"'/^ aeR. (2.5) 

Proof. It is easily verified that M is invariant and its characteristic function satisfies 
(2.1), where F is given by (2.5). □ 
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Corollary 2.7. Let uj denote the Gaussian measure on C with density given by 

TT"^ exp(— I^P), z E C. Given y G M, let M be the image of the Gaussian product 
measure a;®°° on under the following Borel mapping C°° H: 

C^3C^yi-l + eO-BeH. (2.6) 

{Here ^ = (^1,^2, •••) G is regarded as a row vector, so that the {i, j)-entry 
of the matrix B is equal to y{—l + dCj), h j — ■ ■ ■ ■) Then M e and the 
corresponding function F & Ti is given by 

F(a) = , aeM. (2.7) 

1 — lya 

Proof. The invariance of M is obvious. A direct calculation shows that the char- 
acteristic function of M satisfies (2.1) with F given by (2.7). □ 

We shall call the measures defined in Corollaries 2.5-2.7 the elementary ergodic 
measures. 

Note that we could omit the term —1 in (2.6) and then the numerator in (2.7) 
would disappear. However, due to this term, the mean of M is equal to zero and 
the function (2.7) has the property 

p-iya o 

no) = -, — = l--aV + 0(y') asy^O, (2.8) 

1 — lya I 

i.e., the term of degree 1 in y vanishes. This is important for the following con- 
struction. 

Proposition 2.8. The class T\ contains all functions of the form 

^,.,..(a) = e™'/^n^, (2.9) 

where 71 e E and 72 ^ are arbitrary constants, and x = {xi,X2, • • • ) is a sequence 
of real numbers such that "^x^. < 00. 

Proof. Suppose first the sequence x is finite. Then (2.9) is a finite product of 
functions that belong to the class J-'i due to Corollaries 2.5, 2.6, and 2.7. By 
Corollary 2.2, their product belongs to J-'i also. Finally, if the sequence x is infinite, 
then the assumption ^a;^ < 00 and the estimate (2.8) imply that the product in 
(2.9) is convergent for any a G M; moreover, the result is a continuous function. 
Then Corollary 2.4 implies that this function belongs to .Fi. □ 

Comments. 1) The order of the x^s is unessential, so that x = (xk) is rather a 
point configuration (or multiset) than a sequence. 

2) If ^ Ixfel < 00 then (2.9) may be rewritten as 

F,,,,,M = e'^^"""^'^'/'n T^rTTT' 7i := 71 - E^'^- (2-10) 

^ i iXf^a, ^ 

3) The function (2.9) admits a holomorphic continuation to the horizontal strip 
{2; e C I I Im^l < e}, where = sup \xk\. 

4) The parameters 71, 72, x are uniquely determined by the function F^^^y^^^;. 

5) The characteristic function (2.2) corresponding to the function (2.9) can be 
written as 

f{A) = e'-^' ^-^^ ^'-(^ det TT . , AeH{oo). (2.11) 

V 1 — ixkAJ 
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Theorem 2.9 (Classification Theorem). The class is exhausted by the functions 
of the form (2.9). Thus, the characteristic functions of the ergodic U{oo) -invariant 
Borel probability measures on H are just the functions of the form (2.11), where 
7i, 72, a^i, 3^2, • • • CLre real parameters such that 72 ^ and X] ^1 < 

This resuh gives a description of ergodic measures M & M.: any such M is a 
convolution of the elementary ergodic measures constructed in Corollaries 2.5-2.7. 
The proof will be given in §4: we shall derive Theorem 2.9 from a more general 
result, Theorem 4.1. 

Note that the elementary ergodic measures of Corollaries 2.5 and 2.6 are infinitely 
divisible with respect to convolution, whereas those of Corollary 2.7 are not. 

Remark 2.10. The construction of Corollary 2.7 can be generalized as follows. Fix 
k = 1,2,..., consider the space C^^°° of all complex matrices S with k rows and 
infinitely many columns, and equip it with the Gaussian product measure 

^fcxoo _^®oo^...^^®oo_ ^2.12) 

Let z^xi^ . . . ^Xk be any real numbers, and let X = diag(a;i, . . . , x^) denote the 
diagonal matrix of order k with diagonal entries xi, . . . ,Xk- Consider the mapping 

C''''°° 3E^ z-l + E*XE = B e H. (2.13) 

Then the image of the Gaussian measure (2.12) under the mapping (2.13) is an 
ergodic measure Mz-x-i_,...,xk- Its parameters are 

7i = 2; + a;iH \- Xk, 72 = 0, x = {xi, . . . ,Xk,0,0, . . .) . (2.14) 

Mz;xi,...,xki k — 1,2,..., form a weakly dense subset of M. Further, let H^k 
denote the closed subspace of H formed by matrices of rank ^ k. One can prove 
that the measures Mo-xi....,Xk just those measures M G that are supported 
by H^k- Finally, note that the measures Mo.xi,...,Xk with xi = ■ ■ • = Xk are infinite- 
dimensional analogs of the well-known Wishart distributions (see e.g., Muirhead 
[Mu]). 

Remark 2.11. Let H^{n) C H{n) and C H denote the subsets of nonnegative 
definite matrices. Then Hj^ coincides with \im .Hj^(n) and is a closed cone in H. 
One can prove that a measure M £ is supported by if and only if its 
parameters satisfy the conditions 

72 = 0, ^ 0, 0^2 ^ 0, ^a;A;^7i<oo. (2.17) 

Remark 2.12. Let M E M. and let F G be the corresponding function (2.9). 
We may regard (if, M) as a probability space and the matrix elements Bij of a 
matrix B e H as random variables. Let /j, denote the distribution of the real 
random variable Bn; then F is the characteristic function of We know that M 
is completely determined by /U. Let us describe the structure of /U. It follows from 
(2.9) that n is the convolution of a (not necessarily centered) normal distribution 
with a family of distributions possessing characteristic functions of the form (2.7). 
It is easily verified that if j/ > 0, then the distribution with characteristic function 
(2.7) is supported by the half-line t ^ —y and has density 



t^y-^e-y t^ -y. 
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(2.18) 



This is the shifted exponential distribution with variance xp' and mean 0. If the 
parameter y is negative, it suffices to replace t by —t. Thus, /i is the convolution 
of a normal distribution with a family of modified exponential distributions. 

Remark 2.13. Let M E Ai he an ergodic measure for which not all the parameters 
Xk vanish. Then the characteristic function f{A), see (2.11), cannot be factorized 
into a product of factors each of which depends on a single matrix element Aij 
only. This means that M is not a product measure, so that the matrix elements 
Bpq, as random variables defined on (if, M), are not independent on the whole (in 
the particular case of the measures Mz-xi,...,XkJ this is seen from their construction 
via the mapping (2.13)). However, certain matrix elements are independent. For 
instance, it is evident that the diagonal elements are independent. More generally, 
for any ipi,qi), ■ ■ ■ ,{pm,qm) such that i ^ j implies pi ^ pj, qi 7^ qj, pi 7^ qj, 
the matrix elements Bp^q^, . . . , Bp^q^ are independent. This claim can be deduced 
from (2.2). 

According to Proposition 1.5, Theorem 2.9 also gives a complete description of 
spherical representations of the group G{oo). We shall now discuss the possibility of 
extending spherical representations to some topological completions of the inductive 
group G(oo). 

Let us regard U{oo) as a group of unitary operators in the complex coordinate 
Hilbert space £2, and let U{oo) D U{oo) stand for the group of all unitary operators 
in £2- Note that U{oo) is a topological group with respect to the weak operator 
topology (which coincides on unitary operators with the strong operator topology). 
Further, let H{oo)i D H{oo) (respectively, H{oo)2 D H{oo)) denote the space of 
the trace class (respectively, Hilbert-Schmidt) Hermitian operators in £2, equipped 
with the topology defined by the trace norm || • ||i (respectively, by the Hilbert- 
Schmidt norm || • II2). 

One can check that the actions 



where {u,A) 1— > uAu ^, are continuous (cf. Shale [Sha]). Thus we may form the 
semidirect products 



which are topological groups with respect to the corresponding product topologies. 
Note that the inductive limit group G{oo) = U{oo) K H{oo) is contained as a dense 
subgroup both in G(oo)i and 6^(00)2. 

Corollary 2.14. (i) Any spherical representation of the group G{oo) admits an 
extension to a continuous representation of the topological group G(oo)i. 

(ii) It can be continued to 6^(00)2 if and only if the parameter 71 in (2.11) 
vanishes. 

Proof, (i) Indeed, for any function F of the form (2.9), the corresponding function 
/ on if (00), which is defined by (2.11), can be extended to a continuous function 
on H{oo)i. Then the latter function can be extended to a continuous U{oo)-hi- 
invariant function on the group ^(oo)!, and our claim follows. 

(ii) We argue as in (i), using the fact that / is continuous with respect to the 
Hilbert-Schmidt norm if and only if 71 = 0. □ 



U{oo) X H{oo)2 H{oo)2, 



G(oo)i = U{oo) K H{oo)i, 



G(oo)2 = U{oo) K H{oo)2, 
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§3. Approximation of ergodic measures by orbital measures 

Let X he a separable metric space and C{X) be the Banach space of bounded 
continuous functions on X. By a measure on X we shall always mean a Borel 
probability measure. Recall that a sequence {z^n} of measures on X is said to be 
weakly convergent to a measure v (notation: =r- y) if (V') ^n) (V') z^) as n — >• oo 
for any ifj G C{X), see, e.g., [Bi, Pal]. 

Proposition 3.1. There exists a countable set \1/ C C{X) of functions with the 
following property: a sequence of measures z^i, z^2, • • • weakly converges to a measure 
V provided convergence occurs for any function if) . 

Proof. See, e.g., [Pal, Chapter II, Theorem 6.6]. □ 

Let /C(l) C /C(2) C ... be an ascending chain of compact groups, and let /C = 
linj/C(n) be the corresponding inductive limit topological group. We assume there 
is a jointly continuous action (tt, x) ^ u ■ x of the group K, on the space X. Let 
m„ denote the normalized Haar measure on )C{n). Given a point x & X, let m„(a;) 
denote the image of m„ under the mapping u ^ u ■ x, where u ranges over /C(n), 
i.e., mn{x) is the unique AT (n) -invariant probability measure supported by the orbit 
}C{n) ■ X. 

Theorem 3.2 (Vershik [V, Theorem 1]). Let v he an ergodic JC-invariant Borel 
probability measure on X . Then there exists a point x & X such that mn{x) =^ v as 
n — > oo. Moreover, the set of all points x & X with this property is of full measure 
with respect to v. 

Proof. Given if} e C{X), put 

ipn{x) = {ip,mn{x)) = / i){u- x)mn{du), xeX, (3.1) 

JlC{n) 

and let i/j denote the constant function 

i^{x) = (i/;,!^) -1= (^J iPix) i^idx)^ ■ 1 . (3.2) 

To prove the theorem, it is enough to check that ipn ^ fp almost everywhere for 
any i(; e C{X). Indeed, it will follow that, given an arbitrary countable family of 
functions C C{X), there exists a subset X' of full measure such that 

('(/', mn{x)) — > {ip, v) for any •0 G \E' and any x G X' . (3-3) 

Taking as a subset from Proposition 3.1, we shall obtain runix) ^ v for all 
x^X'. _ 

Now we proceed to check that, given '0 G C(A'), we have '^n ^ almost 
everywhere. 

Remark that the V'n are bounded Borel functions. Indeed, given x ^ X and 
n, the function u i— > il){u ■ x) is continuous on )C{n), so that the integral in (3.1) 
can be approximated by its Riemannian integral sums. It follows that ifj^ is a 
pointwise limit of continuous functions, so it is a Borel function. Its boundedness 
is immediate. 
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Now it is enough to verify the following two claims: 

Claim 1. •0^ — •0 in the metric of L^(A', z/). 

Claim 2. For almost all points x & X, the following limit exists 

i^oo{x) := lim Vn(a^)- 

n—^oo 

Indeed, since the ipnS are uniformly bounded, Claim 2 will imply that ipn ~^ 
•i/'oo in L^iXjv), whence t/^oo = V' almost everywhere, so that '4>n ^ '4' almost 
everywhere. 

To verify Claim 1, consider the natural unitary representation of the group K, in 
the Hilbert space L?{X,v). Let Pn denote the orthoprojection in this space onto 
the subspace of /C(n)-invariant vectors, n = 1, 2, . . . . Then strongly converges, 
as n — s> oo, to P, the orthoprojection onto the subspace of /C-invariant vectors. In 
particular, PnV' ~^ m the metric of LP'{X, v). Since the measure v is ergodic, 
the only /C-invariant vectors are the constants (this is the only place when we use 
the assumption that v is ergodic). Therefore, P?/' is a constant function, which 
clearly equals i}). Finally, P^V' = V'n, so that i/'n — m L'^{X, v). 

Now we shall prove Claim 2 by a method similar to the one used in the proof 
of Birkhoff-Khinchine's individual ergodic theorem (see, e.g., Parthasarathy [Pa2, 
§49]). 

Without loss of generality one may assume that ij) is real-valued. For N = 
1,2,..., put 

En ^ En^iI^)^ \x e X\ sup Vn(a;) > o|, (3.4) 

oo 

E^ = E^{'4;) = ixeX\ sup Vn(a;) > o| = U £;iv(V'). (3.5) 

Note that each Ejs[ is a Borel subset, whence E^o also is a Borel subset. 
The following claim is an analog of the maximal ergodic theorem. 
Claim 3. We have 

/ tl){x)v{dx)'^Q. (3.6) 
Indeed, since {En) is a monotone family of sets, it is enough to check that 

I il}{x)v{dx)^^, AT = 1,2,.... (3.7) 

We can write each En as a disjoint union of Borel subsets. 

En = EiN U • • • U Enn, (3.8) 

where 

EmN = {xeX \ i/jmix) > 0, i/jiix) ^0 for m+l^i^N}. (3.9) 
Then it is enough to show that 

/ 'ip{x)u{dx) ^0, l^m^N. (3.10) 
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To do this, we remark that E^n is invariant relative to the action of /C(m). Since 
v is an invariant measure, it follows that 

/ ip{x)i'{dx) = / ip{u ■ x)'u{dx), u & K,{m) . (3-11) 

mN 

Since the function (u, x) '^{u-x) is continuous on /C(m) x EmN-, we may integrate 
the right-hand side of (3.11) over /C(m) (with respect to the Haar measure) and 
then interchange the integrals over K,{m) and over EmN- This yields 

/ i){x)v{dx)= I tl}m{x)v{dx). (3.12) 

By definition (3.9), is positive on EmN, so that (3.12) is nonnegative. 
Thus, we have checked Claim 3. 
Further, for arbitrary real a < b, put 

^ab = {x & X \ lim '4>n{x) < a <b < lim'0„(a;)}. (3.13) 

This is a Borel subset. Let us establish the double inequality 

aiyiXab)> I i^{x)v{dx)^bv{Xab), (3.14) 

which will imply v{Xab) = 0. 

Indeed, Xab is a /C-invariant Borel subset of X. Let us replace X by Xah and 
apply Claim 3 to the functions 

i^' ■.= {i^-b)\x^, and V":=(a-V)U„.- (3-15) 
By definition (3.13) of Xab-, 

E^i,/;') = E^i,/;") = Xab, (3.16) 

whence 

/ ip\x)u{dx) ^0, / 'ilj"{x)u{dx) ^ 0, (3.17) 

which is equivalent to (3.14). 

Finally, applying (3.14) to various couples of rational numbers a < 6, we see that 

lim il^njx) = liiiiV'n(2^) (3.18) 

almost everywhere. This completes the proof of Claim 2 and of the theorem. □ 

We shall apply Theorem 3.2 to X = H and /C = U{oo). Note that the assump- 
tions of Theorem 3.2 are satisfied. Indeed, H is homeomorphic to a separable metric 
space, because it is essentially a copy of M°°, and, further, the action U {oo)xH — > H 
is jointly continuous. 

For n = 1, 2, ... , denote by A4n the set of C/(n)-invariant probability measures 
that are supported by the C/(n)-orbits in the space H{n). These measures will 
be called orbital measures. Since H{n) is contained in H, we may view orbital 
measures as measures on the space H. 
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Theorem 3.3. For any ergodic measure M e M., there exists a sequence {M„ £ 
•M-n} of orbital measures such that =^ M as n — > oo. 

Proof. We shall write uiniB) instead of ninix); here 5 is a matrix from H. By the 
first claim of Theorem 3.2, there exists B E H such that mn{B) =^ M as n ^ oo. 
Consider the projections On'- H ^ H{n) defined at the beginning of §1 and remark 
that 9n{mn{B)) coincides with the orbital measure in H{n) corresponding to the 
matrix 9n{B). Let us take this orbital measure as M„. If /c ^ n, then mn{B) and 
Mn have the same image under dk] fixing k and letting n ^ oo, we see that the 
measures 9k{Mn) on H{k) weakly converge to the measure 9k{M). Now identify 
H with ]R°° and recall that on ]R°° weak convergence of probability measures is 
equivalent to weak convergence of their finite-dimensional projections (see, e.g., [Bi, 
chapter 1, §3]). Applying this to our sequence {Mn}-, we see that M„ =^ M. □ 

It is convenient to analyze weak convergence of measures in terms of charac- 
teristic functions. Assume 1^,1^1,^2,... are Borel probability measures on H and 

/i: /25 • • • denote their characteristic functions. We shall need the following sim- 
ple claim, which again is essentially a well-known fact about the space R°° . 

Proposition 3.4. Weak convergence ^ ^ on H is equivalent to uniform con- 
vergence fn^fon compact subsets in H{oo). □ 

(Note that any compact subset in H{oo) is always contained in H{n) for suffi- 
ciently large n.) 

Proof. We again pass to finite-dimensional projections and then use the fact that 
weak convergence of probability measures on M."^ is equivalent to uniform conver- 
gence, on compact sets, of their characteristic functions. □ 

To obtain Theorem 3.3, we could use, instead of Theorem 3.2, another general 
result, where we have to specialize Q = G{oo), K. — U{oo): 

Theorem 3.5 (Olshanski [03, Theorem 2.5] and [05, Theorem 22.10]). Let G = 
\in^Q{n) be an inductive limit of separable locally compact groups and K, be a sub- 
group of Q. Let )C{n) — K-H Q{n), so that K, = lin^/C(n). 

Then any extreme fC-bi-invariant continuous positive definite function f on Q, 
normalized at unity, can be approached, uniformly on compact sets of the group 
Q, by a sequence {/n} of extreme IC{n)-bi-invariant continuous positive definite 
normalized functions on the subgroups Q{n). □ 

Note, however, that Theorem 3.2, due to its second claim, provides us with more 
detailed information on the approximation process than Theorem 3.5. 

Remark 3.6. Note that Claim 2 in the proof of Theorem 3.2 can also be deduced 
from Doob's theorem on convergence of (reversed) martingales (see [D, Chapter 
VII, Theorem 4.2]). Indeed, denote by the a-algebra of all /C(n)-invariant Borel 
subsets of the space X. We have Bi ^ B2 1^ . . . , so that, by Doob's theorem, 
as n — > 00, the conditional expectation E{il}\Bn) of the bounded Borel function 
ip converges almost everywhere, with respect to z/, to a function, which is the 
conditional expectation £'('0|i3oo), where B^o = f]Bn. To derive Claim 2, we only 
need to show that £{1/; \ Bn) — V'n almost everywhere with respect to z^, n = 1, 2, . . . . 
By definition of conditional expectation, this means that ipn is >Bn-measurable and 
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the following condition holds: 

AeBn ^ / ij{x)u{dx)^ / i;ri{x)iy{dx). (3.19) 

J A J A 

Since V'n is a /C(n)-invariant Borel function, it is B„-measurable. Further, since 
the action of K. on X is jointly continuous, the function (-u, x) ^ '^{u- x) is contin- 
uous on /C X A", hence is a Borel function on /C(n) x A. By Fubini's theorem. 



/ (/ 

JK.{n) \ Ja 



ip{u ■ x)y{dx) \mn{du) = I I / ip{u ■ x)mn{du) \y{dx) . (3.20) 

lK.{n) \ JA J J A \JK{n) J 

The left-hand side of (3.20) is equal to ^ ^■^{x)v{dx), because A and v are /C(n)- 
invariant, whereas the right-hand side is equal to J^il^nix)i'{dx) by the definition 
of ipn- 

§4. Main Theorem 

We shall deal with a sequence {M^ G A4n}, n — 1,2,..., of orbital measures. 
By fn we denote the characteristic function of Mn, recall that 

fn{A)= f e^*-^(^^)M„(dS), AeH{n), (4.1) 

J fir, 

where Qn stands for the U (n)-orbit that carries Let A(n) = (Ai(n), . . . , Xn{n)) 
be the common spectrum of all the matrices B e ^In- Then fin may be spec- 
ified as the orbit containing the diagonal matrix diagA(r2,) with diagonal entries 
(Ai(?i), . . . , Xn{n)). The eigenvalues Ai(n), . . . , A„(n) may be arranged in any or- 
der; it will be convenient for us to separate the positive and the negative eigenvalues 
and to regard A(n) as a double sequence formed by positive and negative eigenval- 
ues, respectively, written in decreasing order of their absolute values: 

A(n) = (A'(n),A"(n)), (4.2) 

where 

A'in) = {X[{n)^X'^{n)^---^0), 

A"{n) = {X'lin) ^ X'^{n) ^ ■ ■ ■ <: 0) . ^ ' 

The possible zero values may be included either in A'(n) or in A"(n); in fact, we 
prefer to view both A'(n) and A"{n) as infinite sequences with a finite number of 
nonzero terms. 

Theorem 4.1 (Main Theorem). Let {M„ e Mn} be an infinite sequence of orbital 
measures defined by a sequence {Qn C H{n)} of U in) -orbits. For n = 1, 2, . . . , pick 

a matrix from the orbit and write the collection K{n) of its eigenvalues as a 
double sequence (4.2). {Note that A(n) does not depend on the choice of B^.) 
(i) Suppose that the following limits exist: 

x> li^^llM^O, 4'=lim^^^0, k=l,2,..., (4.4) 

n— >oo n n^oo fi 



71 = lim - y;(A'fc(n) + A'fc'(n)) = lim - tr (4.5) 

n—^oo ri n—>-oo ri 

k 

^2= lini ;j2 5]((AUn))2 + (A','(n))2)= lim ^tr{B'J. (4.6) 
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Then the measures M„ weakly converge to an ergodic measure M E A4 with the 
multiplicative characteristic function f defined by 



f{A)= n ^(«)' AeH{oo), (4.7) 

aeSpec(yl) 

where Spec (A) stands for the collection {ai, a2, . . . , 0, 0, . . . } of the eigenvalues of 
A and 

\ 1 — IX 1 — zx^a 

here 71 is given by (4.5), the parameters a;^ ^ and x'l ^ are given by (4.4), anc?, 
finally, 

72 = 72-E((4)' + (4'f), (4.9) 

where 72 «s ^fz-yen 6y (4.6). 

(ii) Conversely, if the measures weakly converge to a probability measure on 
H, then the limits (4.4)-(4.6) do exist. 

Comment. It follows from the definition of the parameters x'j^, x'^, and 72 that 
E((4)' + (4)') ^ 72, so that 72 > 0. 

Corollary 4.2. Let {M^ G A^n} be a sequence of orbital measures that weakly 
converges to a Borel probability measure M on H . Then M E Ai. □ 

Derivation of Theorem 2.9 from Theorem 4.1. Let Fe.Fi and let M e Al be the 
corresponding ergodic measure. We must show that F is of the form (2.9). By 
Theorem 3.3, there exists a sequence {M^ G M.n} that weakly converges to M. 
Next, by claim (ii) of Theorem 4.1, the limits (4.4)-(4.6) exist. Finally, by claim (i) 
of Theorem 4.1, F is of the form (4.8) that coincides with (2.9) up to a reordering 
of the points in a; = (a;i, a;2, . . . ) only; we recall (see Comment 1 to Proposition 2.8) 
that we may order the parameters xi, X2, ... in (2.9) in any way. □ 

Outline of the proof of Theorem 4.1. In §5, we establish a preliminary result - we 
expand the orbital integral (4.1) into a series of Schur polynomials. The proof of 
the theorem is given in §6; it is divided into three steps. 
In Step 1, we check that under assumptions (4.4)-(4.6), 

/n(diag(a,0,0, ...)) R, (4.10) 

where diag(- ■ ■ ) stands for a diagonal matrix and F^^^^^^x is given by (4.8) or, that 
is the same, by (2.9). 

In Step 2, we generalize this to arbitrary diagonal matrices: 

k 

/n(diag(ai,...,afe,0,0, ...)) ^ ^^I'^s.^K)' (ai,...,afe) G (4.11) 

This proves claim (i). 

Finally, in Step 3, using a simple trick, we show that conditions (4.4)-(4.6) are 
indeed necessary (claim (ii)). 
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Let us emphasize that to prove Theorem 2.9 only, one could avoid Step 2. How- 
ever, we need this step to characterize the convergent sequences {M^}. 

Remark 4.3. To check (4.10) or, more generally, (4.11), we consider the Taylor 
series decomposition at zero for the left-hand side and show that its coefficients 
tend, as n — > oo, to the corresponding Taylor coefficients for the right-hand side. 
Note that these coefficients are nothing but the moments of the measures. 

Moreover, it follows from claim (ii) and the proof of (i) that whenever a sequence 
{Mn e A^n} weakly converges to a probability measure M, the moments of M„ 
must tend to the moments of M . Thus, in our situation, weak convergence M 
is always controlled by moments — a fact that is not at all evident a priori. 

Such a "moment method" also works in allied classification problems, related to 
characters of U (oo) (see Vershik-Kerov [VK2]) and spherical functions of GL(oo, C) 
(see Nessonov [N]). We conjecture that it can be used for all families of classical 
symmetric spaces. (About spherical functions on infinite-dimensional symmetric 
spaces, see Olshanskii [05] and Pickrell [Pil].) 

Remark 4.4. Assume that in the spectrum A(?i) = (Ai(n), . . . , Xnin)) there are at 
most k nonzero eigenvalues, where k does not depend on n. Then we can prove 
claim (i) of Theorem 4.1 directly, i.e., without using moments. Indeed, let 

lini = a;.^ l^i^n, Ai(n) = 0, i>k. (4.12) 

n—^oo n 

Then the measure M„ can be viewed as the image of the normalized Haar measure 
of the group U (n) under the mapping 

U(n) 9 'u^'fx*diag(Ai(n),...,Afe(n),0,...,0)w = 5 e H{n) C H. (4.13) 

The matrix B can be rewritten as follows: 

S=(S(n))*X(n)S(n), (4.14) 

where S(n) denotes the k x n matrix formed by the first k rows of the matrix 
u & U (n) multiplied by the scalar ^/rl, and 

X(n) = diag(M!),...,M!)). (4.15) 

Now let n — > oo. Then X{ri) X := diag(a;i, . . . , x^), because of (4.12). Further, 
fix m = 1,2,... and regard the k x m matrix formed by the first m columns of 
S(?i) as a random matrix variable (with respect to the Haar measure of U{n)). It 
is well-known that the limit distribution of this matrix is given by the Gaussian 
product measure u'^^'^ (product of km copies of oj, cf. (2.12)), where u stands for 
the Gaussian measure on C specified in Corollary 2.7. This fact is proved, e.g., in 
Olshanski [05, Lemma 5.3]. It follows that the measures weakly converge, as 
n — s> cxD, to the measure MQ-xi,...,xk ^ defined in Remark 2.10. 

Remark 4.5. In general, for a convergent sequence of orbital measures, the eigen- 
values in the spectrum A(n) must grow linearly in n as n — > oo. But if the limiting 
ergodic measure is Gaussian (that is, x = 0), then the order of growth of the 
eigenvalues becomes equal to a/u. Example: let 

A(n) = i^/rri, • • • , Vt^, -Vl^i • • • > -\/7^)> 

" v ' " V ' 

[n/2] [(n+l)/2] 

then the limiting measure is the Gaussian measure Mo^yfi. 
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§5. Expanding spherical functions into series of Schur polynomials 

In this section, we fix n = 1,2, Let (ai, . . . , a„) e (Ai, . . . , A^) e R", 

and suppose 

^ = diag(ai, . . . , On), A = diag (Ai, . . . , An) (5.1) 

arc the corresponding diagonal matrices. We fix A and deal with the t/(n)-orbit 
O C H{n) passing through A. Let /a stand for the characteristic function of the 
invariant probability measure supported by Q. Then 

/a(A)= / e'''^^''^''~'Uu, (5.2) 

Ju{n) 

where du is the normalized Haar measure on U{n). 

Since /a is an entire function of (ai, . . . , an) £ C"^, it admits an everywhere 
convergent Taylor series expansion. But since /A(ai, • ■ • ,cin) is also symmetric in 
ai, . . . , a„, it is more convenient to rewrite this Taylor series as a series of Schur 
polynomials s^, 

/a(^) = /A(ai, . . . , On) = ^ c^s^(ai, • • • , On), (5.3) 

where fj, ranges over the set of all Young diagrams with at most n rows. (About 
Schur polynomials, see, e.g., Macdonald [M].) 

We shall use some standard notation concerning Young diagrams: /i h m means 
that /U is a partition of m (i.e., m equals |/i|, the number of boxes in fi), i{p.) is 
the number of (nonzero) rows in /x, {p, q) E /J, denotes the box of lying on the 
intersection of pth row and qth column, dim fx is the dimension of the irreducible 
representation of the symmetric group S{m), m = that corresponds to the 
diagram fx. 

Theorem 5.1. The coefficients in (5.3) are given by the following formula: 

S = TT — s,^{iXi,...,iXn)=i^^^ TT — s^(Ai,...,An). 

t t ^ + 1-P , ^t'^ n + q-p 

(5.4) 

Proof. Step 1. Let tt be an irreducible representation of U{n), DimTT be its dimen- 
sion, and X be its normalized character: 

, , tr7r(5f) , 
Dim TT 



It is well known that x satisfies the following functional equation: 

/ xiguhu'^) du ^ x{9)x{h), g,heU{n). 

Ju(n) 



(5.5) 



Now let us take as tt the irreducible polynomial representation with highest 
weight 11= (ill,..., /In)- Then 
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where, given a n x n matrix Spec(5r) stands for the collection (2:1, ... , 2;^) of its 
eigenvalues, and Dim„ denotes the dimension of the representation tt of the group 

U{n). 

Substituting (5.6) into (5.5), we obtain 

/ s^{^Y>^c{guhu~'^))du^ s^(Spec(5f)) s^(Spec(/i)) . (5.7) 

Ju{n) -Uim^/i 

By analytic continuation, this formula holds for arbitrary n x n complex matrices 
g and h. 

Step 2. Let us come back to the orbital integral (5.2), where A and A are given 
by (5.1). We may write 

tr(AnA«-i) ^ ^ ^ p^(Spec (iAwA^-^)), (5.8) 



m>0 



where pi{zi, . . . , z^) = + • • • + 2;^ is the first power sum. 

Recall the well-known identity (see, e.g., Macdonald [M, chapter I, (7.8)]): 

Pr(^i>--->^n) = 5Z dim fj, ■ Sfj,{zi,...,Zn). (5.9) 

Using it, we may rewrite (5.8) as follows: 

gUr(AnA«-)^^ ^ ^ (Spec (a«A«-i )) . (5.10) 

Integrating both sides of (5.10) over u E U{n) and applying the functional equation 
(5.7) with g = A, h = iA, we obtain 

(5.11) 

Step 3. For a box (p, q) G /U, let /i(p, q) denote the corresponding hook length. 
Recall the well-known formulas 

dim/x = m! TT / , Dim„//= TT ^^/^ ^ , 
{p,q)&tJ- ip,q)€n 

see, e.g., [M, Chapter I, §3, Example 4]. It follows 

dim n Yj 1 



n 



ml Dim„ n n + q — p 

Substituting this into (5.11), we obtain (5.4). □ 

(After work on this paper was completed, we learned that the argument presented 
above was used much earlier by James, see his survey paper [J, (60)].) 
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Corollary 5.2. The orbital integral (5.2) is given by the following explicit formula: 

fA{A)= e^''^(^-^- Uu = ^—^ --i (5.12) 

Ju(n) V{ai,...,an)V{iXi,...,iXn) 

Here V{- ■ ■) is the Vandermonde determinant, 

V{zi,...,Zn) = Yl i^j-^k), 

and the coordinates ai, . . . , an, as well as Ai, . . . , A^, are assumed to be pairwise 
distinct. 

Proof. The statement of Theorem 5.1 may be rewritten as follows: 

/A(A) = (n-l)!---0! y ^M(Qi,---,Qn)v(iAi,.. zA.) 
^ ^ ^ ^ V (m+^-l)!(/^2 + n-2)!---//,! 

Using the determinant formula for Schur polynomials, we obtain 

V(ai, . . . , an)V{iXi, . . . ,i\n) . . 

(n-l)!---0! -^^^^^ 

^ det [gf ■ det [itXjr^+^-%,^, 



det[af'=]-,=i.det[(zA,)-'=l"- 

mi>--->m„^0 



E 



mi\m2\ ■ ■ ■mn\ 



Since the numerator in the latter expression is symmetric in mi , . . . , m„ and 
vanishes if some of these numbers are equal, we may drop the assumption mi > 
• • • > rUn. Then we obtain 

F(ai,...,a„)l/(?;Ai,...,zAn) 



(n- l)!---0! 



fA{A) 



^ j_ y det[a7'']-,^,-det[(iA,)-'']-,^, 
n\ ^ mil 17121 ■■■ 

mi,...,m„^0 
1 1 

H. m\. • • • llbyi. 

mi,...,mn^O o-,T 

where a = (cri, . . . , cr„) and r = (ti,...,t^) range over all the permutations of 

1, . . . , 71. 

Changing the order of the summation, we obtain 

^ ^sgn(c7)sgn(T)e*"'^i^"i • • •e'"''"^"" = det [e*''^^'=]^fc=i, 
n. 

cr,T 

which completes the proof. □ 

Remark 5.3. Formula (5.12) is well known and can be proved in a number of differ- 
ent ways. For instance, it can be obtained from the Gelfand-Naimark calculation 
of spherical functions on GL{n, C) (see Gelfand-Naimark [GN]) by a passage to the 
limit. Another way consists in using the radial parts of invariant differential oper- 
ators. Yet another method can be found in [BGV, Section 7.5] (we are grateful to 
Michel Duflo for the latter reference) . Note also that writing the above calculations 
in reverse order, we can deduce Theorem 5.1 from formula (5.12). 
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Corollary 5.4. Let us substitute A = diag (a, 0, ... , 0) into (5.2). Then 



/A(diag(a,0,...,0))= / e-(-^--'^^^ du ^ Y: < ^-f'-'-f^^ a^, 

Ju{n) n(n + 1) • • • (n + m - 1) 

(5.13) 

where km is the mth complete symmetric function. 

Proof. Indeed, note that s^(a, 0, . . . , 0) vanishes unless fx = (m), where m = 
0,1,.... Thus the summation is really taken over the diagrams = (m) only. 
For these diagrams, reduces to km and the product g)e/i(^ + 9 ~ P) turns 
into n{n + 1) • • • (n + m — 1), so that we obtain (5.13). □ 

Remark 5.5. The expansion (5.3) of theorem 5.1 may be written in the following 
form, which emphasizes the symmetry between A and A, 

fA{A)= J2 ^''^'f n )s^{X^,...,Xn)s^{ai,...,an). (5.14) 

The symmetry can already be seen from (5.2). 

§6. Proof of theorem 4.1 
We start with the following simple claim. ^ 

Proposition 6.1. Let /, /i, /2, • • • be analytic functions on satisfying the fol- 
lowing conditions: 

(i) /i, /2, • • • are positive definite and normalized at the origin; 

(ii) the Taylor coefficients of f^ at the origin tend, as n ^ oo, to the corre- 
sponding Taylor coefficients of /; 

(iii) The Taylor decomposition of fn converges absolutely and uniformly on n in 
a neighborhood of the origin that does not depend on n. 

Then fn^f uniformly on compact subsets ofM.^. 

Proof. This is a standard exercise. It is clear that fn f i^n Si neighborhood of 
the origin. Consider the probability measures z^i, ^2, . . . on M'^ that correspond to 
the functions /i, /2, . . . by Bochner's theorem. By Paul Levy's classical continuity 
theorem (see, e.g., Shiryaev's textbook [Shi, chapter III, §3]), the measures Un 
weakly converge to a probability measure u. Let / stand for the characteristic 
function of u. Then fn ^ f uniformly on compact sets of M^. Further, (ii) and 
(iii) imply that / is analytic in a neighborhood of the origin, and this implies that 
/ is analytic on the whole space R'^. Therefore, f = f and our claim follows. □ 

We proceed to prove the theorem. The notation of §4 is maintained. 
Step 1. Let us fix a sequence {A(n)} such that the limits (4.4), (4.5), and (4.6) 
exist. Let us abbreviate 

fn{a) = /n(diag (a, 0, 0, . . . )), a e R. (6.1) 



^See also Remark 6.2 at the end of this section. 
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Note that /n(a) is the characteristic function of the probabihty measure '■= 
the image of the measure M„ under the projection 9i : H(n) — > H(l) = M. 
The purpose of this step is to prove that 

hm /n(a) = F^i,^2,a;(a), a e R, (6.2) 

n^oo 

uniformly on bounded sets in R. (Recall that 71 is given by (4.5), 72 is given by 
(4.6) and (4.9), and x = {x'i,X2, . . .;x'i,X2, . . .), where the latter parameters are 
defined by (4.4). The function F^^ x is given by (4.8) or, that is the same, by 
(2.9).) 

To do this, let us expand both the sides of (6.2) into Taylor series: 



fn{a) = J2 (6-3) 

F,,,,U^) = E cl^^a^. (6.4) 

By Proposition 6.1, it is enough to verify the following two claims: First, 

limc(:)=c^), m = 0,1,2,.... (6.5) 

n— >oo 

Second, the series (6.3) converges absolutely and uniformly on n in a sufficiently 
small neighborhood of the origin, i.e., 

Ic^^l^C^C^, m = 0,1,2,..., (6.6) 

where the constants Ci > 0, C2 > do not depend on n. 
By (5.13), we have 

(n) _ hra{iXi{n),...,iXn{n)) 
~ n(n + l)---(n + m-l) 

hJ'-^,...,'^). (6.7) 



n(n + 1) • • • (n + m — 1) \ n n 
It is clear that in both claims, (6.5) and (6.6), we may replace Cm^ by 

(n) ._h f ^-^l(^) i^jn) 



. (6.8) 

Further, instead of the series ^Y^c>m cJ^ and X^c^'^a"^, it is more convenient to 
deal with their logarithms ln(^c^-'a"^) and ln(^ c^''a"^), respectively. 
Now recall a well-known identity from the theory of symmetric functions. 



J] /t^(.)a"' = exp { J]^'^^*)")' 



(6.9) 



where prn( • ) are the power sum symmetric functions (see, e.g., [M]). It follows from 
(6.8) and (6.9) that 
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On the other hand, by definition (2.9) of the function F, 



7i,72,a;) 



In ( J] c(^)a- ) = Z7ia - \ 72a^ + J] p^{ix) ^, 



m^2 



(6.11) 



where 



(6.12) 



k=l 



Note that the sum in the right-hand side of (6.12) is convergent, because, due to 
assumption (4.6), we have P2{x) < oo (see the Comment to Theorem 4.1). 
Thus, our first claim reduces to the existence of the fimits 



'iAi(n) iXn{n) 

hm pi , . . . , 

n^oo \ n n 



iXijn) iXnjn) 
, . . . , 

n n 

iXijn) iXnjn) 

■)•••■) 

n n 



lim p2 

n— »oo 

hm p^ 



-72 + P2{ix), 



= Pmiix), m ^ 3, 



(6.13) 
(6.14) 
(6.15) 



and our second claim reduces to an estimate of the form 

'iAi(n) iXn{n)' 



Pr 



n 



n 



^C;(C^)"^, m^3. 



(6.16) 



where C[ > 0, C2 > are some constants not depending on n. 

Clearly, (6.13) is just the assumption (4.5). Further, (6.14) immediately follows 
from (4.6) and (4.9). Indeed, by (4.6), 



[i\i{n) iXn{n) 
hm p2 I 



n 



hm p2 , . . . , = -72 ■ 

n^oo \ n n 



Now, by (4.9), 



-72 = -72 - J]((4)' + (40') = -12+P2{ix). 



Let us verify (6.15). Using the notation (4.2), (4.3), we have 

'iKin) 



' iXijn) iXnjn) 

Pm \ 1 • • • 1 

n n 



n 



n 



Pmiix) = + E«)' 



(6.17) 
(6.18) 



To deduce (6.15) from (4.4), it suffices to show that both sums in the right-hand 
side of (6.17) converge absolutely and uniformly on n. Let us examine the first 
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sum (for the second one the reasoning is just the same). Since m ^ 3 and A^(n) ^ 
A2(n) ^ • • • , we have for AT = 1, 2, . . . 



i\'^{n) 



n 



r>N 



Kin) 



n 



A'^(n) /Ai(n) 



n 



n 



KM 

n 

An(n) 
n 



m-2 



E 



Kin) 



n 



(6.19) 



By (4.6), p2iXiin)/n, . . . , Xnin)/n) remains bounded as n — oo. Finally, since 
^'n — >■ as — >• oo and since A^(n)/n — >• x'^ as n — >• oo, the value of A^(n)/n may 
be made arbitrarily small provided first AT and then n are chosen large enough. 
Let us verify (6.16). For m ^ 3 



Pr 



iAi(n) iKin) 



n 



n 



^ P2 
^P2 



Ai(n) Xnin) 

■)•••■> 

n n 

Ai(n) Xnin) 



sup 

l<r<r 



Xrin) 



n 



m-2 



n 



n 



P2 



Ai(n) Xnin) 



n 



n 



(m-2)/2 



(6.20) 



Since p2(Ai(n)/n, . . . , A„(n)/n) remains bounded as n — > oo, we obtain (6.16). 
This completes Step 1. 

Step 2. The purpose of this step is to prove that under the same assumptions as 
in Step 1, we have a more general result: for any fixed A; = 1, 2, ... , 



lim /n(diag(ai,...,afe,0, ...,0)) = TT i^7i,72,a;(op) 

p=l 



(6.21) 



uniformly on bounded subsets in M'^. 
Let us abbreviate 



/n(ai, . . . , ak) = /n(diag(ai, . . . , a^, 0, . . . , 0)) 



(6.22) 



and note that /n(ai, ■ ■ ■ ,cik) is again the characteristic function of some probability 
measure i/n on M*^. Namely, i^n is the radial part of the measure 
on H{k) with respect to the projection 



Hik) 3 Spec (A) e 



(6.23) 



Our reasoning will be similar to that of Step 1. We expand both sides of (6.21) 
into multidimensional Taylor series. However, since these are symmetric functions 
of (ai, . . . , ajt), we prefer to rewrite the Taylor series as series of Schur polynomials 
s^(ai, . . . , Ofc), where // ranges over the set of all Young diagrams with £( • ) ^ A;: 



/n(ai, . . . , afe) = ^ cjf )s;^(ai, . . . , a^), 

k 

n ^7i,72,a;K) = c^^hf^iai, ...,ak). 

p=l 
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(6.24) 
(6.25) 



As in Step 1, applying Proposition 6.1, we reduce our problem to verifying the 
following two claims: 
First, 

lim c^') = 4°°) for any /x with ^ k, (6.26) 

n— >oo ^ ^ 

and, second, 

|c(r)KCiCl'^l, (6.27) 

where Ci > 0, C2 > are some constants not depending on n. 
By (5.4), 

(p,<j)€/x 

n f iXiin) i\n{n) 

n -\- 

Thus, in both claims, (6.26) and (6.27), we may replace c}f^ by 



TT ^^-..P^,...,^^ . (6.28) 



Let us prove (6.26) with replaced by ^. 

Recall the Jacobi-Trudi identity expressing the Schur functions in terms of the 
complete symmetric functions (see, e.g., [M, Chapter I, (3.4)]): 

= det[V_i+,]f,^=i, e{fi) ^ k. (6.30) 

It follows from (6.8), (6.29), and (6.30) that 

^=<^\s^^+j\^J=l■ (6.31) 

By Step 1, as n ^ 00, so that 

}i^J^^=^^<7-^^A=l■ (6-32) 

Further, it is well known that for an arbitrary formal series X]^>o Cm.O'^ with cq = 1, 
we have 

n(l]c^<)= 1] det[c^._i+jl^^^-^is^(ai,...,afc). (6.33) 

Applying (6.33) to the left-hand side of (6.25), we conclude that 

= det[c;,7i^,.]t., = ^lim c(r). (6.34) 

Thus, we have verified the first claim. As for the second one, an estimate of 
type (6.27) for c>^^ follows at once from (6.31) and from the estimate (6.6) for Cm"*, 
proved in Step 1. 
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This completes Step 2. 

It follows that fn f (where / is given by (4.7)) uniformly on compact sub- 
sets of if(oo). Clearly, / is a continuous positive definite normalized function on 
H{oo). By Proposition 1.1, it is the characteristic function of a C/(oo) -invariant 
Borel probability measure M on the space H, invariant under the action of U{oo). 
Then, by Proposition 3.4, M„ weakly converges to M. Since / is multiplicative, M 
is ergodic. Thus, we have verified claim (i) of Theorem 4.1. 

Step 3. Let us fix a sequence {A(?i)}, where A(?i) = (Ai(n), . . . , Xn{n)), and let 
fnid), a e M, be defined as in Step 1. Let us assume that 

lim U{a) = f{a), a e R, (6.35) 

n— >oo 

uniformly on bounded subsets in M, where f{a) is a function on M (of course, / 
is automatically continuous). We shall prove that then {A(n)} must satisfy the 
assumptions (4.4)-(4.6) of Theorem 4.1(i). 
Suppose first that 

r /Ai(n) An(n)\ / /Ai(n) A.(n)\y\ 

Then, given an infinite subset A'^ C {1, 2, . . . }, there exists a possibly smaller infinite 
subset N' O N such that the limits (4.4), (4.5), and (4.6) exist provided n goes to 
infinity inside N' . Then, by Step 1, — > F^^^^^^x-, as n ^ oo inside N' , so that 
-^71 ,72, a; — /• For any other N and N' , the parameters 71, 72, and x will be the 
same, because they are uniquely determined by the function itself, see Comment 4 
after Proposition 2.8. It follows that the limits (4.4)-(4.6) exist as n ranges over 
the set of all natural numbers. 

Suppose now that (6.36) does not hold. We shall show that this leads to a 
contradiction with the initial assumption (6.35). 

Indeed, since the expression {• • • } in (6.36) is a homogeneous function of A(n), 
we can choose an infinite subset A?" C {1, 2, ... } and a sequence of positive numbers 
{£„ I n e N} such that lim„gjv s^n = and 

/ / Ai(n) A„(n)\ / / Ai(n) K{n)\V\ 
Inn <^P2 £n ,---,en + Pi ,---,en ( = (6-37) 



neN y \ n n J \ \ n , , , 

Then, replacing A^ by a smaller infinite subset A^', we can arrange so that for the 
sequence {£„A(n)}, the limits (4.4)-(4.6) will exist provided n goes to infinity inside 
N' . Moreover, at least one of the corresponding parameters 71, 72 will be nonzero. 

Note that the effect of multiplying A(n) by Sn is the same as that of multiplying 
a by en- Thus, by Step 1, 

n(i) — -^71,72,2;) ^ I^- (6.38) 

Since at least one of the parameters 71, 72 of the function F^^^^^^x is nonzero, it 
follows from the definition of this function that it is not equal identically to 1. But 
analytic, the same is true in an arbitrarily small neighborhood of 
the point a = 0. Then, comparing (6.38) with (6.35), we arrive at a contradiction. 
Thus, we have verified claim (ii) of Theorem 4.1. □ 
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Remark 6.2. Note that the estimates (6.6) and (6.16) are not necessary to assert 
the convergence of the functions Indeed, Proposition 6.1 may be replaced by 
the following stronger claim: 

Let /i, /2, . . . and f be smooth positive definite functions on 9} , normalized at 
the origin. Expand them into Taylor series at the origin and assume that each 
Taylor coefficient of fn tends, as n —> oo, to the corresponding coefficient of f. 
Finally, assume that the moment problem defined by the coefficients of f has a 
unique solution {the latter condition is satisfied, e.g., if f is analytic). 

Then the sequence {fn) converges to f uniformly on compact subsets o/M"". 

By virtue of this claim, the verification of the uniform convergence of the Taylor 
expansions may be omitted. 

§7. Total positivity 

Definition 7.1. Let (f{t) be a real nonnegative measurable function on M. Then 
(fi is said to be a totally positive function if for n = 1, 2, . . . 

det [(p{ti — Sj)]2j^i ^ for any ti < ■ ■ ■ < tn and si < • • • < . (7.1) 

It will be convenient for us to include in the definition the following additional 
assumption: (p is summable and J (f{t) dt = 1, i.e., (f{t) dt is a probability measure 
on M. (In [S2], functions satisfying both conditions are called Polya frequency 
functions, the second condition is in fact not restrictive, see [S2, lemma 4].) 

Proposition 7.2 (Schoenberg [S2, p. 341, Lemma 5]). The set of totally positive 
functions is stable under convolution. 

Proof. For two summable functions </? and if), the convolution if * ip is correctly 
defined and the following formula is readily verified: 

det[{v * mi -Sj)]lj=^ 

= — det[if{ti-Uk)]lk=i-^^ibP{uk-Sj)]kj^idui---dun. (7.2) 

Now suppose (f and i/j are totally positive and let ti <•■■< tn, si <■■• < Sn- 
Then the integrand in (7.2) is nonnegative for all (pairwise distinct) ui,...,Un, 
because both determinants have the same sign, equal to that of Ylkyii'^k —ui). □ 

Proposition 7.3 (Schoenberg [S2, p. 335 and p. 343]). The densities of the normal 
and exponential distributions, 

V^,(t) = ^e-*'/(^^), 7>0, (7.3) 

and _^ 

are totally positive. □ 

Note that the result remains true after the shift t ^-^ t + const of the argument 
or the change of sign t ^ —t. 

Note also that {'07) and {ipy} are one-parametric semigroups with respect to 
the convolution product. 
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Theorem 7.4 (Schoenberg's theorem on totally positive functions, see Schoen- 
berg [S2], Karlin [K]). The Fourier transforms (p of totally positive functions 99 are 
just the functions F^i,'y2,x, defined in (2.9), where at least one of the parameters 
72,Xi,X2,... is nonzero. □ 

This fundamental result shows that a totally positive function 99 is a convolution 
product ipo * ipi * ip2 * ■ ■ ■ , where v'o(^) dt is a normal distribution and ipk{t) dt, 
k = 1,2,..., are, up to transformations 1 1— > ±t + const, exponential distributions. 

(We have to exclude 72 = xi = X2 = • • • = 0, because the inverse Fourier 
transform of the corresponding function F is a Dirac measure.) 

Comparing Theorem 2.9 and Theorem 7.4, we obtain the following correspon- 
dence M ^ ip between the ergodic measures M e (except the Dirac measures 
on scalar matrices) and the totally positive functions on the real line: 

0i{M){dt) = (p{t)dt 

(recall that the mapping 9i assigns to a matrix B ^ H its matrix element Bn G M). 
In other words, this means that the distribution of the random variable Bn with 
respect to the probability distribution M on the matrices B e H is given by the 

density (f. 

The easy part of Theorem 7.4 consists in verifying the fact that Fy-^^^^^^ = (p with 
a totally positive This is done by making use of Proposition 7.2 and Proposition 
7.3 and an evident passage to the limit (cf. Proposition 2.8). 

The hard part of Theorem 7.4 is to prove that the Fourier transform (p of any 
totally positive function is of the form (2.9). Our purpose is to show that this 
claim is equivalent to Theorem 2.9. 

Definition 7.5 (Karlin [K, p. 49]). A real smooth nonnegative function </7(t) on 
M, J Lp{t) = 1, is called extended totally positive if 

det [cp^'-^^ ivj)]l^^, >0, n = l,2,...,vi>--->vn. (7.5) 

Proposition 7.6. (i) Any smooth totally positive function (p is extended totally 
positive. 

(ii) Conversely, if (p is an extended totally positive function, then the function 
ip * i/jy, where i/j^ was defined by (7.3), is totally positive for any 7 > 0. 

Proof, (i) By definition of total positivity, for any pairwise distinct real ti, . . . ,tn 
and any Si < • • • < 

Uitp - t,)-' ■ det Mt, - s,)]l^^, ^ 0. (7.6) 

p>q 

Putting si = — fi, . . . , Sn = —Vn and letting ti, . . . , ^ in (7.6), we obtain (7.5). 

(ii) By Theorem 2.1 in Karlin [K, p. 50], if <p verifies strict inequalities in (7.5), 
then it is totally positive. So it suffices to prove that if we replace cp hy (p * ip.y , then 
the inequalities in (7.5) become strict. 

For n = 1 this is evident, because (p is nonnegative and not identically equal 
to zero whereas i/j^ is strictly positive. For n > 1 this argument is generalized as 
follows. 
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First, remark that the functions (/?, (/?', . . . are hnearly independent. Indeed, 
if this is not true, then ip satisfies a hnear differential equation with constant co- 
efficients, whence ^ooast— >ooort— > — oo. But this contradicts the 
assumption Lp e L^(M). 

Next, substitute V' = V'l iiito formula (7.2) and repeat the argument used in the 
proof of (i). Then we obtain, for any si < • • • < Sni 

det[(^*V'.)^^-'H-^.-)].Vi 

= - \ det[(^^^-i)(-«fe)]^fe^i det[V'7(«fc - Sj)W^^^ dui ■■■dun. (7.7) 

Arguing just as in the proof of Proposition 7.2, we see that the integrand in (7.7) 
is nonnegative. 

Finally, as noted in Schoenberg [S2, p. 336], the second determinant in the 
integrand is nonzero provided ui, . . . ,Un are pairwise distinct. On the other hand, 
since (p, (fi', ■ ■ ■ are linearly independent, the first determinant in the integrand does 
not vanish for certain (wi, . . . , tt^) G R", hence on an open subset of M". We 
conclude that the integrand is everywhere nonnegative and strictly positive on an 
open subset, so that the integral (7.7) is strictly positive. □ 

Note that Proposition 7.6 corresponds to a part of Pickrell's proof in [Pi2, p. 154- 
155] . There it is claimed that an analytic extended totally positive function is totally 
positive; however, the arguments are too sketchy and seem to be incomplete. To 
avoid this difficulty, we modified the claim somewhat and used a trick suggested by 
Boyer's paper [Bo, p. 218]. 

Note also that not all totally positive functions are smooth: for instance, the 
function (7.4) is not smooth at 0. Thus, the class of extended totally positive 
functions, as defined above, does not coincide with the class of totally positive 
functions (although the two classes are very close, as is seen from Proposition 
7.6). For this reason, to use property (7.5), we must first smooth totally positive 
functions. 

The next theorem is Pickrell's main calculation in [Pi2, pp. 154-155]. It is simple 
but instructive. For completeness and for reader's convenience, we give the proof 
(which is presented here in slightly more detail than in [Pi2]). 

Theorem 7.7 (Pickrell [Pi2, pp. 154-155]). Let ip he a smooth nonnegative function 
on R, / (p{t) dt = 1, and let F = (p he its Fourier transform. Then F helongs to the 
class !F\ {see definition (2.2)) if and only if ip is extended totally positive. 

Proof. Given n = 1,2,... and A e H{n), denote by ai, . . . , the eigenvalues of 
A and put 

/„(A) = F(ai)---F(a^); 

this is a continuous [/(n) -invariant function on H{n). Let us fix n and show that 
positive definiteness of is equivalent to condition (7.5). 

By Bochner's theorem, is positive definite if and only if its inverse Fourier 
transform is a measure. This condition is equivalent to the following one: for any 
function ^ from the Schwartz space S{H{n)), 

{fn, *) := / fn{A)HA)dA ^ 0. (7.8) 

JH{n) 
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Since fn is C/(n)-invariant, one may assume ^ is t/(n) -invariant too. 

Let D{n) denote the subspace of diagonal matrices in H{n). We identify D{n) 
with and write elements of D{n) as diag(ai, . . . , a„) where (ai, . . . , a^) G M"". 
It is well known and easily verified that the radial part of the Lebesgue measure on 
H{n) with respect to the action of U{n) is the measure 

const V^(ai, a„) dai ••• cia„, const > (7.9) 

on D{n), where 

V{ai, . . . , a„) = (ap - a^) . 

p<q 

It follows that 



(/n,$) = const / V{ai,...,an)(p{ai)---(p{an) 

JD{n) 



X F(ai,...,On)*(diag(ai,...,an))(iai---dan. (7.10) 

Put 

e{tu ...,tn)^ V{ti, ^'(diag(ti, . . . , tn)) ■ (7.11) 
We shall show that the integral (7.10) is equal, up to a positive factor, to 

/ det [v^(^-^) {tk)]lk=ie{ti, ...,tn)dti--- dtn . (7.12) 

Indeed, the function V{ai, . . . , a^) <^(«i) • ■ ■ ^(^n) is the Fourier transform of 

V J^, . . . , ■ Vih) ■ ■ ■ V>{tn) = det {tk)]l,^, . (7.13) 

On the other hand, using the t/(n)-invariance of we have 
F(ai , . . . , a„) $ (diag (ai , . . . , On) ) 

= y(ai,...,an) / gitr(diag(ai,...,a„)T)^^y)^y 

JTeHin) 

= V{ai,...,an) [ ([ e^*'^^'^^^s(ai,...,a„)uTn-i)^ 

JHin) \Ju{n) J (7-14) 

Using formula (5.12) for the interior integral and again applying formula (7.8) for 
the radial part of the Lebesgue measure, we see that (7.14) is equal, up to a positive 
factor, to 

(_^)n(n-l)/2 f det[e^'^^*'=]^,fe=i-^(ti,...,t„)(iti---cZtn, (7.15) 

where 6 was defined in (7.11). Developing the determinant and using the fact 
that 9{ti, . . . ,tn) is antisymmetric with respect to permutations of ti, . . . ,tn, we 
conclude that (7.15) is equal to 

n! (-i)^("-^)/2^(ai, . . . , a„) . (7.16) 
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Now (7.13) and (7.16) imply that the integral (7.10) is equal, up to a positive 
factor, to 

(_l)n(n-l)/2 f aet[ip^^-^\tk)]lk=l-0{h,...,tn)dh---dtn 

= / det[ip^^-^\tk)]lk=i-0{tu...,tn)dti...dtn. (7.17) 
Jr" 

Thus, we have verified (7.12). 

Since the integrand in (7.12) is symmetric with respect to permutations of 
ti, . . . ,tn, we see that (7.12) is equal, up to a positive factor, to 

/ det ^(^-1) itk)]lk=i ■ ^iti, ...,tn)dti... dtn . (7.18) 

Recall that 9 is given by (7.11), where \& is an C/(n)-invariant nonnegative function 
from the Schwartz space, and remark that V{ti, . . . ,tn) > in the domain ti > 
■ ■ ■ > tn- It follows that (7.18) is nonnegative for any such if and only if ip is 
extended totally positive. □ 

Corollary 7.8. Theorem 2.9 {classification of ergodic measures M ^ M. or of 
functions F & Ti) and Schoenberg's Classification Theorem 7.4 can be derived one 
from another. 

Proof. Let us show that Theorem 2.9 implies Theorem 7.4. Let <^ be a totally 
positive function. We must prove that the Fourier transform F = is of the form 
(2.9), where at least one of the parameters 72, xi,X2, ... is nonzero. 

First suppose (p is smooth. Then, by Proposition 7.6 (i), (p is extended totally 
positive. Next, by Theorem 7.7, F e JF^, and, finally, by Theorem 2.9, F is of the 
form (2.9). 

The general case can be reduced to that of a smooth (p as follows. We again 
use the Gaussian totally positive function i/j^ (see (7.3)) to smooth (p. Then, for 
any 7 > 0, we have a smooth (even analytic) totally positive function (p*il)^. This 
implies that F{a)e~"'°' is of the form (2.9) for any 7 > 0, whence F itself is of 
this form. 

It should be added that F cannot be equal to F^^,o,05 because we know that the 
inverse Fourier transform of F is a function and not a Dirac measure. So at least 
one of the parameters 72, Xi, X2, ... is nonzero. 

The inverse implication is verified similarly, by making use of Proposition 7.6 
(ii). □ 

Thus, our proof of Theorem 2.9 leads to a new proof of Schoenberg's Theorem 7.4. 

§8. Totally positive functions as limits of splines 

After reading the preliminary version [OV] of the present paper, Andrei Ok- 
ounkov remarked that the one-dimensional projections of orbital measures coincide 
with the so-called fundamental splines (= S-splines) whose limits were studied in 
an important paper by Curry and Schoenberg [CS]. The purpose of this section is 
to briefiy discuss the relationship between the results of [CS] and our results. 

We start by stating some classical facts used in [CS] . 
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Fix real numbers ti < ■ ■ • < tn, called the knots (n ^ 3). There exists a (unique) 
function Mn-i{t) = Mn-i{t; ti, . . . , tn) on M such that: 

(i) On each open interval determined by adjacent knots, M„_i(t) is a polyno- 
mial of degree n — 2. 

(ii) Mn-i^t) vanishes when t < ti or t > tn- 

(iii) Mn-i{t) has n — 3 continuous derivatives at each knot. 

(iv) jMn-i{t)dt=l. 

In [CS], the function Mn-i{t) is called the fundamental spline (with knots 
ti, . . . ^tn). Another term, used in the modern literature, is B -spline. 

The fundamental spline Mn-i{t) is given by the following explicit formula: 

M„_.(M....,u=(»-i)E 'T''\:*i°!r - 

k=l i.l.i^k\ k ''i/ 

Let (Tn-i denote the standard (n — 1) -dimensional simplex, 

an-i^ {ipi,...,Pn) \ O^Pi,...,Pn^l, pi + ■ ■ ■ + = 1} C , (8.2) 

and let ^ denote the affine functional on the simplex taking values ti, . . . ,tn at 
its vertices. Then Mn-i{t) coincides with the density of the image under ^ of the 
Lebesgue measure on a"n-i, so normalized that the volume of Cn-i is equal to 1. 
This implies, in particular, that M„_i(t) is nonnegative, so that Mn-i{t) dt is a 
probability measure on M. 

Using a passage to the limit, one easily extends the definition of Mn-i{t) to the 
case when some of the knots coincide. 

For further properties of the functions Mn-i{t) and for proofs of the facts men- 
tioned above, see [CS] or, e.g., Babenko's textbook [Ba]. 

Now we are in a position to state the main result of Curry and Schoenberg [CS, 
Theorem 6]: 

Theorem 8.1. Consider the class of probability measures on IR which can be ob- 
tained as weak limits of measures of the form Mn-i{t; ti, . . . , tn) dt, where n ^ oo 
and the knots ti, . . . , t„ depend on n. Then the characteristic functions of the mea- 
sures of this class are exactly those given by formula (4.8). 

We shall briefly describe the method of proof used in [CS]. Given ti, . . . ,tn, put 
Fn{a) = J (l- — ) Mn-i{t;t^,...,tn)dt. (8.3) 

Since 

Fn{a) may be viewed as the 'approximate Fourier transform' of M„_i(t). Its ad- 
vantage with respect to the ordinary Fourier image of Mn-i{t) is that it is given 
by a very simple expression, namely 

K(a) = n(l-^)" • (8-5) 

k=l 
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Now note that Fn{ia)~^ is a polynomial in a with only real zeros, and use a 
well-known theorem, due to Laguerre and Polya, which describes the class of entire 
functions that can be approximated by polynomials with real zeros: up to change 
of a variable a i— > ia, these are exactly the reciprocals to functions of type (4.8), see 
Hirschman-Widder [HW]. 

Note that this result on entire functions also plays an important role in Schoen- 
berg's classification of totally positive functions (Theorem 7.4 above). 

Comparing Theorem 8.1 with Theorem 7.4, we conclude that a probability mea- 
sure on M (distinct from a Dirac measure) can be approximated by a sequence of 
fundamental splines with growing number n of knots if and only if it is given by a 
totally positive density. 

The following fact, remarked by Andrei Okounkov, is crucial for our discussion. 

Proposition 8.2 (A. Yu. Okounkov). Let Ai ^ . . . ^ A„ fee real numbers and 
n ^ 3. Consider the U{n)-orbit in H{n) passing through the diagonal matrix A = 
diag(Ai, . . . , A„) and denote by ji^dB) the corresponding orbital measure. 

Then the image of n under the projection H{n) 9 B i— > Sn e M coincides with 
the fundamental spline Mn-i{t; Ai, . . . , A^) dt. 

Proof. For u e U{n), we have 

n n 

(ttAtt~^)ii = ^uikXk{u~^)ki = ^ \uikfXk. (8.6) 

k=l k=l 

When the matrix u ranges over U (n) , its first row 

Z= {zi,...,Zn) := (wii,...,win) (8.7) 

ranges over the unit sphere S'^'^~^ C C"^, and under the mapping u i— > 2;, the orbital 
measure projects onto the normalized invariant measure on S'^'^~^, which may be 
written as 

d(Rezi)d(Imzi) ■ ■ ■d(ReZn)d{Im.Zn) 

— rf(i,,P+...+|,„p_i) ■ m 

Further, under the mapping 

{Zi, ...,Zn)^ • • • > \Znf) = {Pli ■■■iPn) (8.9) 

of the sphere onto the simplex cr„_i, the measure (8.8) projects onto the measure 

dpi ■ ■ ■ dpn , s 

const — -, 8.10) 

C^(Pl + ---+Pn-l) ^ ' 

which coincides with the normalized Lebesgue measure on the simplex. 

Finally, the right-hand side of (8.6), which may be rewritten as Pk^k, is just 

the value at the point (pi, . . . ,Pn) £ Cn-i of the linear functional ^ : R"' — > R taking 
values Ai, . . . , An at the vertices of the simplex. By a property of the spline function 
Mn-i{t) mentioned above, we conclude that the image of the orbital measure 
under the mapping B i— Bn is equal to Mn-i{t; Ai, . . . , A^) dt. □ 

By virtue of Proposition 8.2, the one-dimensional projections of orbital mea- 
sures admit a very nice analytic interpretation in terms of splines, and Curry- 
Schoenberg's result described in Theorem 8.1 turns out to be almost equivalent to 
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the 'one-dimensional part' of our Theorem 4.1, that is, to the results of Steps 1 
and 3 in §6. The difference in the statements is that Curry and Schoenberg do 
not obtain necessary and sufficient conditions on the knots under which a sequence 
of fundamental splines would be weakly convergent; they are only interested in 
describing the limiting functions. 

Theorems 7.4 and 8.1 together imply that totally positive functions are exactly 
those functions which may be approximated by fundamental splines with a growing 
number of knots. This fact seems to be highly nontrivial, because the fundamental 
splines themselves are not totally positive. We hope that the chain of relations 
traced in the present paper furnishes a certain explanation of this phenomenon. 

Remark 8.3. Recall that the Dirichlet distribution V{9i, . . . ,9n) with parameters 
^1 > 0, . . . , ^„ > is defined as the probability measure on the (n — 1) -dimensional 
simplex (8.2) whose density with respect to the Lebesgue measure is given by 

constp^i"^---p^""\ (8.11) 

see Kingman [Ki, Section 9.1]. Now in Proposition 8.2 let us replace the space 
H(n) of n X n Hermitian matrices by the space of n x n real symmetric (respec- 
tively, quaternion Hermitian) matrices. Then the onc-dimensional projections of 
orbital measures coincide with various one-dimensional projections of the Dirichlet 
distribution V{l/2, . . . , 1/2) (respectively, of the Dirichlet distribution V{2, . . . , 2)). 

More generally, one can consider one-dimensional projections of the Dirichlet 
distribution T>{0,...,9) with arbitrary parameter ^ > 0. If ^ = 1, 2, 3, . . . then 
one-dimensional projections of this distribution are the fundamental splines with 
multiple knots: the multiplicity of each knot is equal to 9. For general ^ > 
there is no such interpretation. However, for any 9, using the 'moment method', 
one can still obtain an analog of Theorem 8.1. The limiting measures will have a 
characteristic function of the following form (cf. (4.8)): 

F(a) = /2 TT ^rr^ aeR, (8.12) 

where the parameters are the same as in (4.8). 

Finally, note that the parameter ol = 9~^ exactly corresponds to the parameter 
that appears in the theory of Jack's symmetric functions (see Stanley [Sta] or the 
2nd edition (1995) of Macdonald's book [M]). In particular, the mth moment of a 
one-dimensional projection of the Dirichlet distribition T>{9^ ■ ■ ■ ,9) is equal, up to 
a scalar factor, to Pm{ti, ... ,tn', 9~^), where ti, . . . ,tn stand for the parameters 
of the projection (i.e., the values of the corresponding affine functional ^ at the 
vertices of the simplex) and Pm{ • ; a), m = 1, 2, . . . , are one-row Jack's symmetric 
functions with parameter a. 
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